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We report on a phenomenological study of superfluid to Mott insulator transitions of bosons on 
the triangular lattice, focusing primarily on the interplay between Mott localization and geometrical 
charge frustration at 1/2-filling. A general dual vortex field theory is developed for arbitrary rational 
filling factors /, based on the appropriate projective symmetry group. At the simple non-frustrated 
density / = 1/3, we uncover an example of a deconfmed quantum critical point very similar to 
that found on the half-filled square lattice. Turning to / = 1/2, the behavior is quite different. 
Here, we find that the low-energy action describing the Mott transition has an emergent nonabelian 
SU(2) x (7(1) symmetry, not present at the microscopic level. This large nonabelian symmetry 
is directly related to the frustration-induced quasi-degeneracy between many charge-ordered states 
not related by microscopic symmetries. Through this "pseudospin" SU(2) symmetry, the charged 
excitations in the insulator close to the Mott transition develop a skyrmion-like character. This 
leads to an understanding of the recently discovered supersolid phase of the triangular lattice XXZ 
model— as a "partially melted" Mott insulator. The latter picture naturally explains a number 
of puzzling numerical observations of the properties of this supersolid. Moreover, we predict that 
the nearby quantum phase transition from this supersolid to the Mott insulator is in the recently- 
discovered non-compact CP 1 critical universality class^- A description of a broad range of other 
Mott and supersolid states, and a diverse set of quantum critical points between them, is also 
provided. 
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I. INTRODUCTION 

Theoretical interest in quantum phase transitions from 
superfluid to Mott insulating states of bosons has re- 
cently been revived by their experimental realization in 
cold atoms trapped in an optical lattices 1 Extensions of 
these experiments are expected to soon provide a great 
variety of real life toy models, where theoretical sce- 
narios for such phase transitions can be testedA^ For 
the condensed matter community, such transitions are 
interesting in their own right, but also provide a sim- 
pler context in which some aspects of Mott conducting- 
insulating transitions of electrons can be explored. A 
better understanding of such Mott criticality generally 
may help to explain mysteries in various strongly corre- 
lated materials, from heavy fermion metals^ to cuprate 
superconductors^ in which Mott criticality may plausi- 
bly be argued to play a key role. 

An exciting theoretical development in the field has 
been the discovery that some quantum phase transitions 
require a fundamentally new description, not based on 
the now-standard Landau's concept of an order param- 
eter. Instead, such quantum critical points (QCPs) are 
described in terms of emergent degrees of freedom, not 
present in either phase and appearing due to certain spe- 
cial dynamically generated low-energy symmetries at the 
critical point. 7 An interesting consequence of the emer- 
gent low-energy symmetry of the critical point is the 
near-degeneracy of "competing ordered" states unrelated 
by any microscopic symmetry (but unified with one an- 
other by the emergent one) in the neighborhood of the 
quantum phase transition. There is, at present, unfortu- 
nately, no general way to a priori identify the appropriate 
emergent degrees of freedom, should they exist, for any 



putative quantum critical point. 

In the particular context of two-dimensional bosonic 
superfluid-insulator transitions, a general non-Landau- 
Ginzburg- Wilson (non-LGW) framework has recently 
been proposed in Ref. |8] (and see Ref. |9| for a pedagogical 
review), and carried out explicitly on the square lattice. 
In particular, the Mott transition can be described in 
terms of the vortex excitations of the superfluid. In a 
two dimensional superfluid, vortices are point-like "par- 
ticles" whose creation/annihilation operators can be used 
to construct a quantum field theory. The vortices being 
non-local topological objects, these vortex fields are, how- 
ever, not themselves order parameters in the LGW sense. 
The non-locality is manifested by the presence of a (non- 
compact) U(l) gauge field, to which the vortex fields are 
coupled in the "dual" vortex field theory. This formula- 
tion is general because it is based on the excitations of the 
superfluid, which is a stable and apparently featureless 
(i.e. without broken symmetry apart from off-diagonal 
long range order) state at any boson density. Never- 
theless, it was shown in Ref. that the vortices exhibit 
a subtle quantum order which is sensitive to the boson 
density / (per unit cell of the lattice). In particular, at 
non-integral /, the vortices form non-trivial multiplets 
transforming under a projective symmetry group (PSG 
- technically, a projective representation of the lattice 
space group). The Lagrange density for the vortex field 
theory therefore takes the general form 

N 

C = ^2[\( d ^- iA M 2 + S \ ( Pi\ 2 ] +Ant[{^}] 



(1) 



where cpi with 



1 ... N are vortex fields for the N 
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members of the multiplet (N depends upon / - see 
Sec. |nj, and is the dual U(l) gauge field. Quartic 
and higher order terms are contained in £i n t, the struc- 
ture of which is dictated by the PSG. The Mott transition 
is captured by Eq. Q in a simple way. The ground state 
of s > corresponds to the vacuum of vortices - flux 
is expelled from the system, so this is a superfluid. In 
the s < phase, at least one of the vortex flavors will 
condense, and the gauge fields will acquire a gap by the 
Higgs mechanism. This is indicative of a charge gap, and 
describes an incompressible Mott insulating phase. 

Transformations within the ip£ multiplet comprise the 
emergent symmetry operations of the "deconfined" quan- 
tum critical point. This structure has an important phys- 
ical consequence as well: it naturally and unavoidably 
leads to (particular) broken spatial symmetries in the 
Mott state. A mean-field analysis of the vortex field the- 
ory predicts a direct superfluid to Mott transition, as well 
as the nature of the charge ordering in the Mott phase. 

In this paper, we extend this vortex field theory ap- 
proach to describe superfluid to Mott insulator transi- 
tions of bosons on the triangular lattice (hexagonal lat- 
tice in proper crystallographic nomenclature) at frac- 
tional boson fillings / = p/q, with p, q relatively prime 
integers. We focus particularly on the most interesting 
examples, / = 1/2 and / = 1/3. The 1/2-filling case 
introduces a new ingredient not present on the square 
lattice: geometrical frustration. Here we refer to "charge 
frustration" of the ordering of localized boson configura- 
tions in the presence of short-range repulsive interactions. 
A consequence - or perhaps definition - of such geomet- 
rical frustration is the exact or near degeneracy of many 
distinct ordered states. The similarity of this property 
with the emergent near-degeneracy of competing orders 
near a deconfined quantum critical point suggests a pos- 
sible link between the two phenomena. This connection 
indeed appears to be borne out by the analysis in this 
paper. 

In the simplest classical models of frustration, the de- 
generacy amongst low-energy states is not only large but 
macroscopic (i.e. with an entropy proportional to the 
sample volume). This classical degeneracy, when lifted 
by quantum fluctuations, may produce unusual ground 
states. A number of very recent papersifiaii 1 ^ have in- 
vestigated the system of hardcore bosons with nearest- 
neighbor interactions (which can be mapped to a spin- 
1/2 XXZ model) on the triangular lattice near half-filling. 
This realizes such approximately classical frustration in 
the limit of very strong near-neighbor repulsion. These 
works demonstrated that in this system the lifting of the 
macroscopic classical degeneracy results in an unusual 
supersolid ground state, which we denote SS3 because of 
its 3-sublattice structure. This "order by disorder" mech- 
anism is very different from the "conventional" (theoret- 
ically!) picture of supersolidity, via a condensation of 
vacancies and/or interstitials in an ordered solid. 13 

Since the above studies clarified that such Mott states 
do not occur in the simplest nearest-neighbor interac- 



tion model, it is apparent that microscopically, longer- 
range interactions, possibly including ring-exchange, 14 
are necessary to observe these transitions in a microscopic 
model. Finding simple interactions that produce nontriv- 
ial insulating ground states on the triangular lattice is an 
important and difficult problem, that will likely require 
sophisticated numerical methods. We will not address 
this issue here (but see the Discussion, Sec. IVjl . 

Our phenomenological vortex field theory, on the con- 
trary, describes universal aspects (independent of micro- 
scopic realization) of Mott insulating, superfluid, and 
other states, and the transitions between them. As for 
the square lattice case, a mean-field analysis predicts a 
direct superfluid-Mott transition, with a diverse set of 
Mott insulating phases. Also like the square lattice, the 
vortex field theory has an enhanced emergent symmetry 
at the critical point, a hallmark of deconfined critical- 
ity. A significant difference from these prior examples of 
deconfined criticality, however, is that, in the frustrated 
case, / = 1/2, the emergent symmetry is nonabelian, 
containing an SU (2) "pseudospin" subgroup. The much 
larger (than in non-frustrated cases) emergent symme- 
try can be understood physically as symptomatic of the 
larger near-degeneracies present in this case due to frus- 
tration. 

Remarkably, going beyond the mean-field analysis of 
the vortex field theory, our approach connects very nicely 
to the supersolid phase of the XXZ model. Indeed, 
the supersolid order parameter - describing the growth 
of the supersolid state out of the featureless superfluid 
- appears in a particularly simple form in the vortex 
variables. Moreover, our approach reveals an alterna- 
tive view of the supersolid, as a partially-melted "par- 
ent" Mott insulating state, with "quantum disordered" 
pseudospin. This loss of pseudospin order simultane- 
ously with the onset of superfluidity is possible because, 
as we show, the pseudospin skyrmion excitation of the 
Mott insulator carries physical boson charge. The su- 
persolid may thereby also be viewed as a condensate of 
these skyrmions. Furthermore, this picture leads directly 
to the prediction that the Mott insulator to SS3 transi- 
tion is described by the recently discovered Non-Compact 
CP 1 (NCCP 1 ) quantum critical universality class^This 
transition describes the quantum disordering of a pseu- 
dospin vector S (the order parameter for the additional 
solid order of the Mott state) in 2 + 1 dimensions, when 
"hedgehog" instantons are absent in space-time. These 
hedgehog events correspond precisely to processes which 
change the skyrmion number, and are therefore prohib- 
ited by charge conservation in this case. 

The paper is organized as follows. In Sec[n]we develop 
the dual vortex theory for the triangular lattice at a gen- 
eral rational filling. We derive the PSG transformations 
of the degenerate low energy vortex modes and discuss 
some of their general properties. In Sec. lIIII we apply this 
theory to the two cases / = 1/3 and / = 1/2, and discuss 
the resulting Mott states that are obtained by a mean 
field analysis of the vortex field theory. In Sec. IIVI we 
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present a "hard-spin" formulation of the vortex action, 
which enables a study of its phases and transitions be- 
yond mean field theory. We discuss the new phases which 
arise, notably supersolids, the elementary excitations of 
the different states, and the transitions amongst them. 
This includes notably the SS3 phases and the NCCP 1 
transitions between them and their parent Mott insula- 
tor. We conclude with a discussion in Sec.^of a more mi- 
croscopic physical picture for the most interesting Mott 
and supersolid states at half-filling, the connection to the 
recent studies of the XXZ model, and the prospects for 
observing more of the physics in this paper in related 
models. A number of appendices provide useful details 
of various technical results. 



II. CONTINUUM DUAL VORTEX THEORY 

As discussed in the introduction, our aim is to derive 
a field theoretic description of the vicinity of a superfluid 
to Mott insulator transition, in terms of the vortex exci- 
tations of the superfluid state. As is well-known, lattice 
models of bosons can indeed be reformulated in vortex 
variables on the dual lattice, a technique called duality. 
As detailed for instance in Ref. 9 (and references therein), 
this can in principle be carried out exactly for any lat- 
tice boson Hamiltonian, e.g. Bose-Hubbard models, XXZ 
models, etc. Unfortunately, this exact mapping does not 
lead to a particularly tractable limit of the vortex the- 
ory, and it is therefore difficult to extract quantitative 
predictions directly from the lattice vortex theory. 

Fortunately, we can avoid this difficulty in addressing 
our stated goal of understanding universal phenomena 
in the vicinity of the superfiuid-Mott quantum critical 
point. For this purpose, we need not specify any par- 
ticular microscopic boson model. Instead, we will illus- 
trate the calculations through the use of a very simple 
lattice vortex theory, which is chosen to have the same 
spatial symmetries as the physical triangular lattice, and 
to exhibit a superfluid to Mott insulation transition. The 
universal properties of interest will coincide with those of 
more physical microscopic boson models. 

In the Euclidean vortex coherent-state path integral 
formulation in which we work, the lattice vortex action 
is 

S = ^^Ke-^i+.lc.c] 

a /.i 

+ $>IV>.| 2 + t*hk| 4 ] 

a 

+ ^^(^AA^ aA -27r/V) 2 - (2) 

Here a labels sites of the dual 2+1-dimensional uniformly 
stacked honeycomb lattice, with the vertical direction be- 
ing the imaginary time direction r, and fi is summed 
over nearest-neighbor links (3 at 120 degrees connect- 
ing spatial neighbors and a fourth along the imaginary 




FIG. 1: Dual honeycomb lattice. Reflection axes di and d,2 
are shown by dashed lines. 

time direction). The action is written in terms of the 
complex vortex field ?p a (and its conjugate ^*) which an- 
nihilates (creates) unit vorticity, as well as a dual gauge 
field A afl . The physical meaning of the gauge field A afi 
is that its curl gives (2tt times) the bosonic 3-current 
density, and importantly, the temporal component of the 
current density is (2tt times) the charge density. For a 
bosonic system with density / bosons per site, we must 
therefore enforce the condition that 2wf flux on average 
passes through each hexagonal plaquette of the honey- 
comb lattice. We will assume the filling to be rational 
/ = p/q> where p and q are relatively prime integers, and 
will mostly concentrate on the two cases / = 1/2,1/3. 
Physically, t v represents a vortex hopping amplitude, s 
and u represent short-range vortex "core" energies and 
interactions, and e 2 represents the strength of the dual 
electromagnetic field fluctuations (it is roughly propor- 
tional to the local superfluid density away from the Mott 
quantum critical point). 

We may think of the parameter s in Eq.© as driv- 
ing the superfluid-insulator transition. Large positive s 
corresponds to the superfluid phase, in which vortices 
are gapped, and large negative s corresponds to the set 
of possible insulating phases, in which vortices are con- 
densed. It is convenient to think about the transition in 
Eq. @ for small e 2 , i.e. neglecting to a first approxi- 
mation dual gauge fluctuations - they will however be 
restored at a later stage of analysis. In this limit we can 
treat the dual gauge field in a mean field approximation 
and take: 

e^AA„A aA = 2tt/(5 mt . (3) 

As s is decreased from large positive values, we expect 
an instability of the "vortex vacuum" when the energy of 
the lowest vortex excitation approaches zero - or equiv- 
alently, when the minimum eigenvalue of the quadratic 
form for ip a , ^* in Eq. J5J vanishes. To study the univer- 
sal critical properties of the superfluid-insulator transi- 
tion, it is sufficient to isolate these low-energy particles, 
which comprise the vortex multiplets discussed in the in- 
troduction. The continuum limit of Eq. @ then consists 
of a set of vortex fields, each representing one of these 
minimum energy vortex particles. We will also take the 



4 



trivial continuum limit of Eq.(0l in the temporal direc- 
tion. 

Because of the non-zero gauge flux through each spa- 
tial plaquette, the minimum energy multiplets are non- 
trivial. The form of the continuum action in this case is 
determined by the projective representation of the space 
group (PSG)jS under which the vortex fields ip transform. 
To find it, we must work through the consequences of 
some specific gauge choice. As in Ref|$, we will choose 
the Landau gauge for A a/1 . Namely, let 



under 



ai 



1 

x, a 2 = --x- 



V3. 



-v, 



(4) 



be the two basis vectors of the honeycomb lattice, as 
shown in Fig^ (note that the honeycomb lattice has two 
sites per unit cell). Coordinates will be specified, when 
explicit, in this basis, r = aiai + a2a2, with integer oi, a 2 , 
and a\ = a 2 = corresponding to a "type 1" site (see 
Fig. ^) of the dual honeycomb lattice. Landau gauge is 



A ay = 2nfai, 



(5) 



and A n 



for all other directions; that is, only the 



gauge field on vertical links is chosen non-zero. The full 
PSG is generated by a set of unitary transformations of 
ip, one for each generator of the lattice space group. For 
the triangular lattice, we choose the space group genera- 
tors as two elementary lattice translations Xi,T 2 , a 27r/3 
rotation with respect to site 1 in Fig. R 27r / 3 , an d two 
reflections, 1^,1^*^ The PSG transformations of the 
vortex fields, characterized by the unimodular complex 
number ui = e 2lTl f are then: 



R 



Ti 

2tt/3 



ipi{ai,a 2 ) -> -0i(ai-l,a 2 )w a2 

V> 2 (ai,a 2 ) -> ip 2 {ai~l,a 2 )uJ ai+1 ' 

i/>i(ai,a 2 ) -> i/ , i(ai,a 2 -l) 

^2(^1,02) -> 4> 2 {ai,a 2 -l) ' 

^1(01,02) -> ^1(02-01,-01)0;"^^ 

^2(01,02) -> ^2(02-01,-01-1 

Maua 2 ) -» Vi*(-a2,-ai)^ ai ° 2 

Mai,a 2 ) -» ^(-a 2 -l,-ai-l)^ ai(Q2+1) 

^i(ai,a 2 ) -> ip 2 (a 2 ,ai-l)uj aia2 

ih(ai,a 2 ) -» ^i*(a 2 + l,aiV ai(a2+1) ' 



l(2o 2 -ai + l) 



C : ^,-(01,02) -> i/ , l *(ai,a 2 ), 



(7) 



and simultaneous sign change of the fluctuating part of 



the dual gauge field A^ 



-A„ 



The quadratic action of Eq. in Landau gauge has a 
periodicity in real space of q unit cells in the ai direction, 
and one unit in the a2 direction (note that, of course, the 
physics itself has the full periodicity of the honeycomb (or 
underlying triangular direct) lattice). The eigenstates of 
the quadratic action can therefore be characterized by 
their quasimomenta in the corresponding reduced Bril- 
louin zone. Specifically, we introduce the basis vectors of 
the reciprocal lattice, 



bi 



V3 



y- 



v/3 



(8) 



Wavevectors will, when necessary, be specified by coor- 
dinates (ki, k 2 ), with k = fcibi + k 2 b 2 (reciprocal lattice 
vectors correspond to k\ , k 2 being integral multiples of 
2tt). 

By applying the methods of Ref|8|, one may readily 
find the minimum energy multiplet for arbitrary p, q, and 
their PSG transformations. Briefly this is accomplished 
by Fourier transforming Eqs. ijfJJl to obtain the PSG for 
general i/>(k) in momentum space - this is given in Ap- 
pendix ^ - and using the non-commutative algebra of 
translations and rotations implicit in Eqs. © to gener- 
ate a full set of eigenfunctions. One finds two cases. For 
q odd, there are q minima of the vortex dispersion, i.e. 
TV = q in Eq. Q). These occur at momenta 



k< = (0,27r/£), e = 0,...,q-l, 



(9) 



taking wavevectors k^ to lie in the reduced Brilloin zone 
—Tt/q < k\ < w/q and — 7r < k 2 < tt. By contrast, for q 
even, N = 2q, and it is convenient to divide minima into 
two sets a = ±1 = ±, parameterizing I = (1 + a)q/2 + a, 
iff = ... q — 1. The vortex field operator (p acr then acts 
on eigenstates with the wavevector 



(ira/3q, ira/3q + 2irfa) 



(10) 



(6)in the reduced Brillouin zone. For even q, therefore, 
Eq. may be rewritten as 



For the special case / = 1/2, in which we are principally 
interested, one usually considers in addition to these spa- 
tial symmetries, an extra particle-hole symmetry, which 
we denote C for "charge conjugation" . This interchanges 
singly-occupied and empty sites on the direct lattice. In 
the spin language appropriate for hard-core bosons, this 
is just a 180° rotation around the x axis in spin space, 
which has the effect of an Ising transformation S* — > —S* 
(and likewise for Sf). The XXZ model, and indeed any 
hard core boson model at / = 1/2 with only pairwise 
interactions, possesses such a particle-hole symmetry. In 
the dual theory, this requires the action to be invariant 



q-l 



£ even = £ £ [| ( ^ _ ^ ^ |2 + ^ |2] 

i{e^xd v A x f. (11) 



a=± cr=0 
1 

2?' 



We note in passing that the set of ipi corresponding to 
these wavevectors gives the smallest number of minima 
possible for the vortex dispersion (i.e. they comprise the 
smallest-dimensional representation of the PSG for the 
given /), and they are what is realized without special 
fine-tuning of the vortex kinetic energy terms. 
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The PSG transformations of the multiplet can be found 
using Ea. HAll) of Appendix 1X1 by a straightforward gen- 
eralization of the procedure, detailed in Ref. ||. For in- 
stance, under the translations, one finds 



for q odd, and 
T 2 



T 2 



Pao 
facr 



-lj 



(12) 



(13) 



for g even. Here and in the following, the index I for odd 
q and a for even q will be regarded as cyclic modulo q. 
The remaing PSG generators are given in Eqs. IA2IA3I 
in Appendix IaI 

It is now straightforward to write down the most gen- 
eral continuum vortex Lagrangian density, describing the 
superfluid-insulator transition on the triangular lattice. 
The most general terms at quadratic order are simply 
given by Eqs. (|llll|) . Let us now consider the quartic 
potential terms in the continuum theory. Following the 
general approach of Ref. 8, we first write down the con- 
tinuum Lagrangian, imposing only the restrictions from 
gauge symmetry and translational symmetry, i.e. invari- 
ance under Eqs. (|12I13|) . One finds 



9-1 



CI 



odd 



^ImnVlVt+mVl+nVt+m-n-i ( 14 ) 
£mn 

9-1 

Ect@ * * 



a.ai ,0-2—0 



and 7°^/ an 



for odd and even q, respectively, where 7. 
arbitrary at this stage. 

Imposing additional restrictions on the coefficients in 
the above Lagrangians from invariancc under rotations 
and reflections and taking into account hermiticity and 
permutation symmetries, one may obtain a set of condi- 
tions on the 7 coefficients required to preserve the full 
triangular lattice symmetry. These conditions are given 
in Eqs. i|A8IA7|) . For specific /, they can readily be solved 
to derive explicit forms for the Lagrangian. We will give 
these explicitly for / = 1/2 and 1/3 below. 

Interestingly, it is possible to make at least one general 
observation concerning the symmetry properties of £f on . 
All terms in the continuum Lagrangian possess a global 
vortex U(l) symmetry, 



9a 



Pa 



(15) 



That is just a consequence of gauge invariance, expressing 
the conservation of the bosonic current. However, it is 
clear from Eq. H14[) that the quartic Lagrangian £| ven 
possesses (at least) another, "staggered" U(l) symmetry: 



This emergent U(l) symmetry of the vortex theory im- 
plies that there are (at least) two conserved dual charges, 
that can be labelled by the index a. As discussed in 
Ref. |S|, this is linked to the appearance of fractionally- 
charged bosonic excitations at the critical point. We will 
elaborate on the nature of these excitations in Sec. II VI 



III. (DUAL) MEAN FIELD THEORY 

In this section we will discuss a mean field analysis of 
the vortex theory, focusing on the nature of the ordered 
Mott insulating states that occur. We first present some 
general aspects of how spatial order parameters are con- 
structed in the vortex formalism, and give some physical 
picture of how to think of the different Mott states. The 
remaining two subsections describe the specific phase dia- 
gram in the case / = 1/3 and the much more complicated 
and more interesting case / = 1/2. 



A. Order parameters and Mott states 

General argumentsiSiii and physical reasoning seem 
to imply that, barring exotic situations such as phases 
with "topological order" ,!& Mott insulating states occur- 
ing for non-integral / must break space group symmetries 
(and in particular translational symmetry). Such space 
group symmetry breaking is measured by spatial order 
parameters, the simplest of which (sufficient for our pur- 
poses) describe non-uniformity of the boson density (be- 
yond that which is imposed by the underlying triangular 
substrate) . 

To visualize ordering patterns in the insulating phases 
we will find it convenient to introduce a general "density" 
function p(r), where r is a continuous real-space coordi- 
nate with r = taken to coincide with a honeycomb 
lattice site of type "1" (Fig. QJ. We construct p(r) to 
have the property that it transforms like a scalar boson 
density under all symmetry operations. It will be con- 
venient to plot p(r) to graphically illustrate the symme- 
try of the non-uniform states that emerge in the theory. 
Writing r = r\a.i + r^s^-, one can actually construct such 
a function quite generally for odd values of q: 



q{t) = 



m.n 



mri +nr2 



(17) 



where the Fourier components g mn serve as order param- 
eters for different ordered states, and are given by: 



9-1 



b{m,n)Lo " ) lo fe^e+n- (18) 



=0 



</W — ► </we 



(16) 



S(m, n) here is a scalar form factor, that can not be deter- 
mined from symmetry considerations, but should be cho- 
sen to depend only upon the magnitude of the wavevec- 
tor mbi + n\>2- A convenient and simple choice is the 
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Lorentzian, 



S(m, n) 



l + m 2 + (m + 2n) 2 /3' 



(19) 



which we use only for plotting purposes. It is easy to 
check that Qmn indeed transform like Fourier components 
of density: 



T\ . Qmn ^ 
-*2 '• Qmn - * 
-R2-7T/3 : Qmn 
Idi ■ Qmn > 

Id,2 ■ Qm 



Qmn i 
Qmn t 
^ Qn, — 7n—n-> 

Q—n, — mi 

(n-m)/3 
W Qnm- 



(20) 



A similar function can be constructed for q even. When 
/ = 1/q, it takes the form 



Q(r) = E [fi 



Q-n 



a 2-Kia(r l +r 2 )/3q 



(21) 



where the density wave amplitudes are given by 



S(m, n ) e 7riQ ("^ m )/ 3, ? t j^ m "/ 2 +(™- m )/6 

9-1 

5>" m V a V?+n. ( 22 ) 



9-1 
£=0 



£=0 



Here the amplitude S(m, n) has been taken in a simple 
form consistent with rotational symmetry: 



S(m, n) 



1 + (m + a/3) 2 + (m + 2n + a) 2 /3 



(23) 



Eqs. (|22|l are written in terms of ^72(0;), which also enters 
the PSG in general for even q, see Eqs. I|A3IA5|) of Ap- 
pendix E] It is hard to determine for general q. However, 
for the case we will focus upon, / = 1/2, one has 



V2(a)=- — 



for / = 1/2. 



(24) 



We will present plots of p(r) for various mean field (and 
beyond, in the following section!) Mott insulating states. 

These images, and their deconstruction into the 
Fourier amplitudes (g m n etc.), characterize the broken 
spatial symmetry of the Mott states. They do not, how- 
ever, directly give a physical picture for the ground state 
itself. Of course, for a general interacting boson model 
(and certainly in our approach where we do not specify 
the microscopic Hamiltonian), we cannot hope to write 
down the ground state wavefunction. Moreover, this has 
far too much information. What is conceptually useful is 
to understand how to write down a simple wavefunction 
appropriate to a Mott insulator with the same symme- 
tries as predicted by the phenomenological theory. By an 



appropriate wavefunction, we mean one explicitly with 
the correct boson filling, with no long-range correlations 
beyond that of the Mott state, and consistent with in- 
compressibility. We consider a satisfactory form to be a 
product state, 



(25) 



The meaning of Eq. I|25[) is as follows. We divide the 
lattice into a set of non-overlapping identical clusters of 
sites - unit cells of the Mott state - labeled by the □. 
The Y[ indicates a direct product over states defined on 
the Hilbert space of each box, with the same state \tfj(0)) 
chosen on each box. The state must be an eigen- 

state of the total boson number on the given cluster, for 
this wavefunction to represent a Mott insulator, and this 
number should be chosen to match the filling, i.e. equal 
to / times the number of sites in the cluster. Clearly this 
state has only local (within a cluster) charge fluctuations, 
consistent with incompressibility. 

Of course in most cases many such wavefunctions can 
be constructed for a Mott state of a given symmetry, and 
moreover the true ground state wavefunction for a real- 
istic hamiltonian will not have the direct product form. 
However, the above general construction can serve the 
purpose of providing an example of a state with the same 
symmetry properties as predicted by the phenomenolog- 
ical theory, and is expected to be in a sense (which we do 
not attempt to define precisely) adiabatically connected 
to all ground states in the Mott phase. To our knowl- 
edge, all well-established examples of bose Mott insulat- 
ing ground states in the theoretical literature (e.g. on 
the square lattice, checkerboard, stripe, and VBS states) 
have such a construction. We therefore view the existence 
of such a wavefunction as a sort of consistency check on 
our results, and give examples in the following subsec- 
tions. When the meaning is obvious, we give only a brief 
physical description of the state, with the understanding 
that one should keep the block product form, Eq. 
in mind. 



B. / = l/3 

At 1/3-filling, the Mott state is not "frustrated" in the 
colloquial sense. The analysis below reveals that this case 
is extremely analogous to the superfluid-Mott transition 
on the square lattice at half-filling, and in particular pro- 
vides another example of a possible deconfined quantum 
critical point of the type discussed in Ref. |7|. 

Solving Ea. (|A7|l at / = 1/3, we obtain the follow- 
ing form of the quartic potential in the continuum La- 
grangian density: 



d =u(|^o| 2 + |vi| 2 + |^ 2 | 
«{^o| 2 |^i| 2 + ^i| 2 |^ 2 | 2 + 



|2|,„_|2 ,, 2Re 



m\ \fo\ +' 



X (<^oVi^2 + <P2 2 t/Wl + <P*l 2l P2<P0 + C.C.)] } 



(26) 
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As discussed in Ref. 0, it is convenient to transform to a 
different set of variables that realize a "permutative rep- 
resentation" of the PSG. As in the square lattice case, 
such representations exist only for special fillings. In par- 
ticular, it is easy to show that a permutative representa- 
tion does not exist at / = 1/2, in contrast to the square 
lattice case. At / = 1/3 however, it does exist and is 
given by: 

Co = -j= (<A) + e- 27r4/ Vi + W>) , 

6 = 77= (e" 2W Vo + <pi + <p 2 ) , 

6 = -^ (<Po + + e- 2 " /3 ^ 2 ) • (27) 

The physical meaning of the three vortex quantum num- 
bers in the permutative representation of the PSG is that 
they represent three conserved (as shown below) dual 
charges (vorticities). The quantum numbers that are 
dual to these three vorticities are three real fractional 
(1/3 of the boson charge) charges. That is, a dual "vor- 
tex" in which the phase of any one of the & fields winds 
by ±27r at infinity carries a localized charge (boson num- 
ber) of ±1/3 (see Ref. for a simple derivation of this 
result in a more general context). 

The representation of the PSG realized by the ii fields 
is permutative in that each symmetry operation is real- 
ized as the composition of a permutation of the ii fields 
and a simple phase rotation. In particular, 

Ti : it — > ie+i, 



Tn it- 










-R27173 ; 


io -> e**'%, ii - 


-> — «6oi 6 






Id! ■ & 




-e" /6 ei, 


6 




Id 2 ■ £0 






6 


(28) 



The quartic potential simplifies greatly in these variables: 

A=^(|£o| 2 + I6I 2 + I6I 2 ) 2 
+ v(\io\ 2 \ii\ 2 + \ii\ 2 \i2\ 2 + \i2\ 2 \io\ 2 ) (29) 

At the quartic level, it is immediately apparent that the 
microscopic overall U(l) gauge symmetry required by the 
vortex non-locality has been elevated to a £7(1) 3 sym- 
metry under independent rotations of all three it fields. 
Of this, only the group of equal rotations of all fields is 
gauge, leaving an additional [/(l) 3 /[/(l) = U(l) x [7(1) 
global (not gauge) symmetry of the dual theory. This 
emergent symmetry is broken at 6 th order by the term 

£2 = w [(Coii) 3 + (£i6) 3 + (do) 3 + c.c] . (30) 

The structure of the theory is thus rather similar to the 
vortex theory at f — 1/2 on the square lattice&SS It 
provides another example of deconfined criticality, if, as 




FIG. 2: Charge density pattern at 1/3-filling for v > 0. 




FIG. 3: Charge density pattern at 1/3-filling for v < and 
w < 0. 



seems likely, the mean-field irrelevance of the higher- 
order term in £ 2 remains valid with fluctuations in 2 + 1 
dimensions, for some sign of v. The mean field phase 
diagram of the vortex theory can now be easily found 
analytically. 

Let us now proceed with the mean field theory for / = 
1/3. For s < and \w\ < \v\ one obtains 3 distinct 
insulating phases: 

I. v > 0: 

The energy is minimized if only one of the 3 vortex fla- 
vors condenses. This state is then clearly 3-fold degen- 
erate and breaks all symmetries except reflection with 
respect to d,2 axis. In Fig|5] this state is visualized ex- 
plicitly by plotting the corresponding density wave or- 
der parameter. This state is the simplest CDW state at 
1/3-filling, with an example wavefunction consisting of 
a boson number eigenstate on each site. Microscopically 
this phase is the natural Mott insulating state in a model 
in which the Mott transition is driven by strong nearest- 
neighbor repulsive interactions, like the XXZ model. We 
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FIG. 4: Charge density pattern at 1/3-filling for v < and 
w > 0. 

do not expect the superfluid-Mott transition to this state 
is likely to be continuous or deconfined when fluctuations 
are taken into account in 2 + 1 dimensions. 

II. v < 0: 

It is energetically favorable to condense all 3 vortex fla- 
vors, so that all vortex fields have equal magnitude. 
There are then 2 different phases, depending on the sign 
of w. These states are those expected to be connected 
by a deconfined quantum critical point to the superfluid 
state. 

1. w < 0. 

Writing ~ |£|e ie<! , the minimum is achieved 
when: 

6»i - O = 27rm/3,6» 2 - 6 = 2im/Z, (31) 

for m,n — 0,1,2. This state is thus 9-fold de- 
generate and corresponds to period-3 site-centered 
stripes, running in ai, &2 and ai +a2 directions, see 
FigGH One may construct a wavefunction for this 
state by e.g. taking linear 3-site clusters at a 60° 
angle to the stripe (e.g. horizontal in Fig. 0, and 
putting a boson at the center of the cluster, leaving 
the other sites empty. 

2. w > 0. 

One can readily verify that the ground state is 
18-fold degenerate, the distinct solutions being ob- 
tained from 

0i - 9o = 9 2 - 9 Q = 2tt/9, (32) 

by applying translations and the reflection I^ 2 . 
The corresponding characteristic density pattern is 
shown in Fig@] It is adiabatically connected to a 
"bubble" phase or crystalline state in which one 
boson is placed on each site of each elementary tri- 
angle (3 bosons total) on a 3 x 3 triangular superlat- 
tice. The 18-fold degeneracy results from 9 states 



obtained by translating this pattern, and another 
9 states obtained by choosing say down-pointing 
instead of up-pointing triangles. 

C. / = 1/2 

1. Action and symmetries 

We now turn to the more complicated and interest- 
ing case of / = 1/2. It is straightforward to show (by 
a simple generalization of the argument used in Ref. 8 
to prove the absence of a permutative representation in 
the / = 1/3 case on the square lattice), that for this 
case there is no permutative representation of the PSG. 
Unlike in the square lattice case, there are also, surpris- 
ingly, more low-energy vortex modes for / = 1/2 case 
than for / = 1/3. The emergent low-energy symmetry 
amongst these modes is, moreover, nonabelian, with an 
SU(2) "pseudo-spin" subgroup that we will uncover be- 
low. This structure has interesting consequences for the 
phases and excitations that will be explored here and in 
the following section in some detail. 

The quartic potential can be found by solving Ea. l|A8|) . 
It turns out to have the following form: 

£i = "fc^) 2 + v{<p* + „<p +<r ){ip*_ a ,<p_„,) (33) 
+w (|(^_ | 2 |v+o| 2 + f-i^+if-of+o + (0 <-> 1)) • 

To make the symmetries of the Lagrangian more trans- 
parent, it is useful to rewrite it as follows. First we define 
z aa pseudospinor variables, 

f+o- = z a , (34) 

where e CT(T < is the antisymmetric tensor with eoi = — eio = 
1. The quartic action has SU(2) symmetry under rota- 
tions of the a index of z aa . This is made manifest by 
introducing the pseudospin vector 

&a = Z aa T aa' Z aa' > (35) 

where r a , a = x, y,z are the Pauli matrices, and sum- 
mation on the repeated a, a' indices is implied. The 
transformation properties of the S± vectors are particu- 
larly simple, and given in Appendix [5] In terms of these 
pseudospin variables the quartic potential becomes, after 
a trivial redefinition of v and w couplings: 

C\ = u (5+ + S-) 2 + v S+S- + w\ S + ■ S_, (36) 

where S a — |S Q | = ^z* aa z aa (sum on a implied). It 
is now clear that the quartic potential has, in addition 
to the microscopic gauge U(l) symmetry, an SU(2) x 
U(l) x Z2 invariance. The SU(2) symmetry is mani- 
fest in Eq.ljSU), the extra U(l) symmetry is the "stag- 
gered" U(l), already mentioned above, see Ea. (jlfjfl . also 
manifest since 5 1 " are independent of the staggered U(l) 
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phase. The Zi symmetry is the interchange S+ <-> S_ 
(particle-hole symmetry C in fact requires this invariance 
up to a sign, though Eq. (j^H is obtained without using 
C). The pseudospin variables are directly related to the 
Q^ nn density components: 

q% x = ae™ {1+a)/6 S*, 

Qii = SI (37) 

The physical import of the SU(2) symmetry of Eq. i|36|) 
is now clear: all CDW states, related to each other by 
arbitrary rotations in the space of the three Fourier com- 
ponents of the density, are degenerate at this order. This 
degeneracy will be weakly broken, of course, by higher 
order terms in the action. As discussed in the introduc- 
tion, the large emergent symmetry is thereby connected 
with geometrical charge frustration at / = 1/2. 




FIG. 5: Ordering wavevectors in the reciprocal lattice. Filled 
symbols are the wavevectors referred to in the text. Open 
symbols indicate wavevectors belonging to the same star 
which are related to the former sets by complex conjugation. 



2. Order parameters 

The pseudospin vectors S± serve as gauge-invariant or- 
der parameters to characterize the breaking of the SU (2) 
symmetry. It is instructive to construct two other such 
order parameters. The emergent U(l) symmetry is best 
characterized by a complex order parameter, tp, defined 
by 



ztt/4 * 

e -\-<J — (T 3 



(38) 



where we have included the 7r/4 phase factor for later 
convenience. One may also define an Ising order paramer 
$, with 



$ = |* + l 2 -k-l 2 



(39) 



Non-zero ($) implies Id 2 and C are broken. The physical 
meaning of non- vanishing S± and "0 will be elucidated in 
detail in the following. 

Though the terms which break these symmetries are 
small near the Mott QCP (hopefully irrelevant there), 
they are important at sufficiently low energy. We must 
therefore consider those higher order terms in the action 
which are required to reduce the SU(2) x U(l) symme- 
try to only what is required by the PSG. Higher order 
terms which do not reduce this symmetry need not be 
considered. 

First consider the SU (2) pseudospin symmetry. In 
the absence of microscopic particle-hole invariance, it is 
broken at 6 th order by a term of the form S+S+S+ + 
S^S^Sl. We will, however, for concreteness focus on 
the case relevant to the XXZ model and other pairwise 
interacting boson lattice models, in which particle-hole 
symmetry is an invariance of the theory. In this case, 
the SU (2) symmetry is broken only at the 8th order by 
a "cubic anisotropy" term, 



£ 2 = ™ 2 £[(^) 4 + (^) 4 + (^) 



(40) 



The "staggered" U(l) symmetry is more persistent and 
is broken at the 12th order. The simplest term at this 
order that breaks the "staggered" U(l) symmetry is: 



C 3 



-w 3 Re (V 6 



(41) 



It is useful to reorganize various terms in the density 
expansion of Eq. (|21|l to understand in more detail the 
nature of the different order parameters. The wavevec- 
tors referred to in the following are labelled in Fig. El 

Consider first Ising order. Non-zero (<fr) implies only 
that the reflection Id 2 and particle-hole symmetry C are 
broken. Thus it corresponds only to a modulation of the 
density within the primitive unit cell of the triangular 
lattice. The corresponding density modulation therefore 
occurs entirely at reciprocal lattice vectors. The smallest 
set of reciprocal lattice wavevectors that can describe the 
modulation are B = (0,27r), B = (2ir, — 2ir), B 2 = 
(— 27r,0). This density modulation is 



p$(r) = $ ^2 cos(B„-r + ^) 



G 



(42) 



While /;$ = on triangular (direct) lattice sites (as it 
must), it alternates sign on sites of the dual honeycomb 
lattice, i.e. takes opposite signs on centers of triangles of 
the direct lattice. 

Now let us turn to pseudospin ordering. The existence 
of non- vanishing (S±) implies charge ordering at the 
wavevectors G x = (0, — 7r), G y — (— tt,tv), G z — (n,0), 
which lie at the centers of the zone edges. In particular, 
the associated density modulations take the form 

p S (r)= J2 Re \(Sle-™/ 3 - S°:y G «- r ] , (43) 



neglecting higher harmonics which do not change the 
symmetry of ps(r). 
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FIG. 6: Charge density pattern at 1/2-filling for v < 0, Wi < 
0, W2 < and W3 < 0. 



Next consider the XY order parameter ip, A state with 
(-0) exhibits a three-sublattice structure, characterized by 
the zone boundary wavevectors Qo = 2w(%,— i), Qi = 



27r(-i|), Q 2 = 27r(-|, -i): 



p^,(r) = Re 



-i7r/6 



e i27rn/3 e iQ»T 

n=0,l,2 



(44) 



In fact, these three wavevectors differ only by reciprocal 
lattice vectors. From this, it is straightforward to show 
that py,(r) vanishes on dual lattice sites, so that all tri- 
angular plaquettes of the direct lattice arc equivalent up 
to rotations in an XY ordered state. 

Finally, simultaneous breaking of pseudospin and XY 
symmetry is characterized by the "composite" order pa- 
rameter d, a complex vector, defined as 



(45) 



When (d) =^ 0, density modulations appear at the 
wavevectors K x = 2ir(±,-±), K y = 27r(-A,-|), K z = 
2ir(— i, i), which lie within the zone. The corresponding 
density is 



p d (r) = Re 



o 37Ti/4 



J2 d Q e 2,r(a - 1)l/ V K - r 



a=0,l,2 



, (46) 



where we have identified a = x,y, z with a = 0, 1, 2 re- 
spectively. 



3. Mean field phases 

The mean field phase diagram of Cq + L\ + £2 + £3 
can be easily obtained analytically. One finds 10 differ- 
ent phases: 2 for v > and 8 for v < 0. We will not 
attempt to be exhaustive in describing these states. We 
will, however, discuss some of the phases in detail, and 



go into some general aspects of the 8 cases with v < 
in Sec. I1VI Minimizing the mean-field energy functional, 
v < implies that both a — ± vortex pseudospinors are 
condensed with equal amplitude, so S+ = 5_. The sign 
of w\ determines the relative pseudospin orientation. For 
w\ < 0, they are parallel, i.e. S + = S_ = S. In terms 
of the vortex variables this condition is most generally 
solved by 



Z±<T 



(47) 



where S — z^Tfjfji z^i . For w\ ^> 0, the two pseudospin 
vectors are antiparallel, S+ = — S_ = S, which implies 



Z+a 
Z-a 



z a e 



i6/2 



,i$/2 



(48) 



We note that the dual gauge symmetry acts differently 
in the two cases. For parallel pseudospins, a gauge trans- 
formation takes Zfj — > e lx z a , while for antiparallel pseu- 
dospins, instead 8^8 + 2%. 

The remaining terms, u>2 and W3, fix the remaining 
non-gauge symmetries. The sign of W2 chooses the easy 
axes of the pseudospin vector, along (100) and symmetry- 
related axes for w 2 < 0, and along (111) and related axes 
for W2 > 0. The above conditions leave only the rela- 
tive phase between the spinors z±, corresponding to the 
"staggered" U(l) symmetry of the corresponding terms 
in the vortex Lagrangian. 

This relative phase is fixed at 12 th order in z± a , for 
instance for w\ < by the £3 term (for wi > 0, a more 
complicated 12 th order term must be included as the W3 
interaction vanishes in that case). Depending on the sign 
of u>3, the energy minimum is achieved when sin(6#) = 
±1, i.e. 



7T 

12 



Tin 



n = 0, . . . , 5, 



n = 0, . . . , 5, 



4 3 

One of the states, corresponding to 8 = 7r/12 and 



(49) 



(50) 



Z +Q — z +i 



Z-0 



Z-l 



71"' 



(51) 



is shown in Fig|H| We will elucidate the physics of this 
particular pair of states in some detail in the discussion 
section. 

We can generally classify all the mean field states by 
their degeneracies and the corresponding unit cell sizes. 
The states with the pseudospin easy axis along (100) 
and parallel pseudospins (the state in Fig|S] belongs to 
this group) are all 36-fold degenerate and have a 6-site 
unit cell. States with parallel pseudospins but with the 
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easy axis along (111) are 48-fold degenerate and have the 
largest, 12-site unit cell. In the case of antiparallcl pseu- 
dospins, ground state degeneracies are the same, but the 
unit cell size of the 36-fold degenerate states doubles to 
12 sites. States with the smallest unit cells are obtained 
when only one of the pseudospinors is condensed. In this 
case one obtains a 6-fold degenerate ground state and a 
2-site unit cell when (100) is the easy axis, and an 8-fold 
degenerate state with a 4-site unit cell when the easy axis 
is along the (111) direction. 

For the phases of most interest (v < 0) in which 
the pseudospin vectors are either parallel or antiparal- 
lel, some of the density functions associated to the order 
parameters in Sec. IIII (T2l can be simplified to an extent. 
The only qualitative case is the pseudospin vector density 
ps- When S + = S_ = S, it reduces to 

p s \\(r)= J2 S a cos(G a -r- T ). (52) 

a—x,y,z 

Most interesting, when S+ = — S_ = S, one has instead 
Psd*)= E S Q cos(G a .r-I[). (53) 

a—x,y : z 

In this case, it is noteworthy that pgw vanishes on all 
triangular lattice sites. This is a consequence of the fact 
that a configuration of antiparallcl pseudospins preserves 
particle-hole symmetry, so modulations can occur only in 
bond or plaquette "kinetic" terms. ThusS + — S_ maybe 
considered a (particular) purely valence bond solid order 
parameter. 

The composite order parameter d also simplifies once 
the pseudospin order is determined. For parallel pseu- 
dospins, one simply has d = S. For antiparallel pseu- 
dospins, instead, one has d x d* = 2iS, which implies 

d = ei — ie-2, e\ x e2 = S, (54) 

so that ei, e2, S form a right-handed orthogonal frame in 
the 0(3) spin space. The angle of ei in the XY plane is 
arbitrary, and determined by (twice) the phase of z a . 

IV. "HARD SPIN" DESCRIPTION: BEYOND 
MEAN-FIELD THEORY 

In the preceding section, we have followed a Landau- 
theory like procedure (albeit with non-LGW vortex 
fields) in expanding the effective action in a power se- 
ries in the ipe fields, whose amplitude is viewed as small 
in the vicinity of the Mott QCP. In low-dimensional sta- 
tistical mechanics, it is often preferable to formulate the 
theory in terms of "hard spin" variables, in which the 
amplitude of the order parameter field(s) is fixed, and 
only the "angular" degrees of freedom are free to fluc- 
tuate and vary in space and time. Examples include 
the Kosterlitz-Thouless theory of the XY phase transi- 
tion, and the non- linear sigma model formulation of 0(n) 



models. The intuitive rationale for such an approach 
is that, in low dimensions, fluctuations suppress the or- 
dering point of the transition well below the mean-field 
point, so that substantial amplitude is already developed 
in the true critical region. Whatever the rationale, there 
are some advantages to such a hard-spin approach. Dual- 
ity transformations generally apply to hard-spin models. 
Hard-spin variables are particularly appropriate to de- 
scribe the elementary excitations of "ordered" phases in 
which the amplitude of the fields is on average large, and 
only the Goldstone-like fluctuations of the orientation of 
these fields comprise low-energy excitations. Finally, a 
hard-spin formulation returned to the lattice is fully- 
regularized, and can thereby address non-perturbative 
phenomena in a controlled manner. 

In this section, we will provide and analyze a hard- 
spin formulation of the dual vortex action. These allow 
us to identify the excitations and their quantum num- 
bers within the Mott phases. Notably, we find that the 
most interesting Mott states support two distinct kinds 
of excitations. First, there are 1/2-charged ("spinon" in 
the spin- 1/2 XXZ language) "vortex" excitations which 
are linearly confined in pairs deep in the Mott state, but 
are only logarithmically interacting up to a long "confine- 
ment length" near the superfluid-Mott transition. Sec- 
ond, there are additional unit charged "skyrmion" ex- 
citations which are everywhere deconfined in the Mott 
state. These are adiabatically connected to single bo- 
son vacancies/interstitials in the Mott solid, but become 
topological as the Mott transition is approached. The 
hard-spin models also provide a firm ground on which 
to study other phases, notably supersolids, which do not 
occur within a mean field treatment of the vortex field 
theory, but are extremely natural in this approach. 



A. Formulation of hard spin model 

To write down an appropriate hard-spin model, we 
imagine tuning s < to a point beyond the mean field 
Mott transition point. At such a point, the minimum ac- 
tion configurations have non-zero amplitude. We will as- 
sume the mean amplitude is determined by a balance be- 
tween the quadratic Lagrangian, Eq. (|ll|l and the quar- 
tic terms in C%, Eq. (|36() . The minimum action saddle 
points of the combination of these two terms are constant 
in space-time, but allow for a continuous set of orienta- 
tions in the field space. We focus in particular on v < 0, 
so that the magnitude S+ = S- is fixed. We will focus 
primarily on the case w± < 0, so that the saddle point 
has S + = S_. This is the most interesting case, because, 
as we shall see, the recently-determined supersolid phase 
of the XXZ model can be understood in this framework. 
At the end of this section, we will briefly summarize the 
results of similar analysis for w\ > 0, corresponding to 
anti-parallel pseudospins. 

We further suppose that s is negative enough that fluc- 
tuations in the above conditions may be neglected, but 
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that within these constraints the z aa fields can vary spa- 
tially. We will therefore absorb any effects of the mag- 
nitude of the fields into coefficients, and without loss of 
generality normalize to S+ = S— = 1. Ultimately, the 
higher order terms will still be included, but can be con- 
sidered small perturbations. 

The most general solution of S + = S_ = S constraint 
in terms of the vortex variables is given by Eq. (|47|l . 
We will therefore rewrite the action in terms of the CP 1 
field z a and an XY field e %d l 2 . Note that this solution 
possesses a Z 2 gauge invariance under 



9 + 2n. 



(55) 



This is in addition to the physical symmetries of the 
model. It is a gauge invariance since it can be performed 
independently at each space-time point without changing 
z Qcr , and hence physical quantities. 

Inserting Eq. l|T7| into the action, and regularizing it 
on a space-time lattice, we obtain 



C 



HS 



1 



with z* a z ia 



-t v e- iAi "z^z i+fiff cos{A f t8 l /2)+^(e f , l/X A„A lX ) 2 , 

(56) 

1 normalized on each site i of the space- 
time lattice. Note that Eq. H56[) is indeed invariant inde- 
pendently under Eq. (|55|l at each point i. It is convenient 
to rewrite the first term in Eq. (|56[) . making the Ising 
gauge symmetry explicit by introducing an Ising gauge 



field Oi^ which resides on the link 



C\ 2 = -t z a ifi e z* ia z i+tia - tga ifl cos(A M 6' l /2) 



+ ^2 ( e ^v\^vAi\) 



(57) 



where we have introduced two independent parameters 
t z ,tg. One could add a plaquette interaction (line prod- 
uct of <7j„ around plaquettes), which is a standard "ki- 
netic term" for the Z 2 gauge field for further generality. 
We are, however, most interested in the limit in which 
it is absent. In this case, one may sum over on each 
link independently, to return to an action of the form 
of Eq. I|56[l (with additional higher-order terms) . We will 
not attempt, however, to constrain Eq. i|57[l to be exactly 
equivalent to Eq. Ij56|l . Instead, since we are anyway con- 
structing a phenomcnological theory, we regard the free- 
dom to vary t z , tg independently as a means of capturing 
the different possible tendencies due to fluctuation effects 
and details of microscopic dynamics in different physical 
systems. 



B. Phase diagram for parallel pseudospins 

Let us now discuss the phase diagram of Eq. (|57|l . For 
t z ,tg <C 1, both 9i and Zi a variables are disordered, and 
the dual gauge field is gapless. This is the superfluid 
phase. For tg ~ t z 3> 1 large and of the same order, we 
expect that both the CP 1 and XY variables are ordered. 




FIG. 7: Possible charge ordering patterns in the SS3 phases. 
Left: "antiferromagnetic" supersolid (11)3 > 0). Right: "ferri- 
magnetic" supersolid (103 < 0). 



This is the Mott insulator, whose precise nature depends 
upon the anisotropy terms we have neglected to write. 

Now suppose tg ^ 1 but i 2 « 1. In this limit, we 
expect the XY variables condense. This is a "Higgs" 
phase 19 from the point of view of the Z% gauge variables: 
the linear coupling to cos(A M #i/2) means that the 
fields can be regarded as having some non-vanishing ex- 
pectation value in this state. The CP 1 fields however 
remain uncondensed, and A^ remains gapless, so this 
state retains superfluidity. It does, however, break spa- 
tial symmetry. To understand the nature of this symme- 
try breaking, let us take for simplicity t z = 0, and again 
imagine "summing out" the Z 2 gauge fields while varying 
tg. The effective Lagrange density is then 



Cxy = V(A M 0i 



^(e^AA.A.A) 2 (58) 



with V(Q) 



In cosh 



V2 



y/l + C0S(8) 



This some- 



what unconventional gradient term has all the same sym- 
metries as the usual cos(A M #i) term in an XY model, 
and indeed reduces to that form for small tg. On in- 
creasing tg, therefore, a 3D=2+1 dimensional XY transi- 
tion is expected, into a state with a non-zero expectation 
value of e . Comparison with Eq. I|38|l indicates that 
if) ^ e ITr / 4 e - * e is an order parameter for this transition. 
Note that e ie ^ 2 is not an order parameter since it is not 
(Z2) gauge invariant. As noted earl ier, ij) is exactly the 
order parameter identified in Refs. lRllllTll2l as charac- 
terizing the supersolid phase of the XXZ model on the 
triangular lattice. This supersolid has a three-sublatticc 
structure, so we will denote it by SS3. Actually there 
are two different ordering patterns possible within the 
tripled unit cell, depending upon the sign of the 6-fold 
anisotropy term, W3, which should be added to Eq. 158|l . 
They are shown in Fig. [7] 

Finally, consider similarly the situation when tg ^ 1 
but t z varies from small to large. For large t z , we then 
expect the CP 1 variables order but the XY variables re- 
main uncondensed. Analogously to the previous case, 
imagine increasing t z from small to large with tg = 0. 
One can again integrate out the Ising gauge field to ob- 
tain an effective action which is i?2-gauge invariant. In 
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this case, there are two distinct types of "kinetic" terms 
which arise on nearest-neighbor bonds. For small t z , they 
take the form 

£su(2) = — *sSj • Sj - t 2 e 2lA ^z* a z* icr ,z i+ll c ,z i+li (T , 
+ ^(e (U/A A 1 ,4 lA ) 2 , (59) 

where ts^2 ~ t 2 z and Sj = z*TZ t . Two distinct types 
of "orderings" are clearly possible on increasing t z . The 
most natural possibility, driven by ts, is for S to order. 
As above, this is a supersolid, but with a different set 
of possible charge order patterns, characterized by the 
zone boundary center wavevectors rather than those at 
zone corners. An alternative possibility, driven by £2 5 
is that the vortex pair field Zi a Zi a i condenses. Such a 
paired vortex condensate is a spin-liquid insulator (see 
Ref. I2H) (a "Z2" spin liquid in the now-conventional 
nomenclatur e) 22 ! 23 ). The gauge fluctuations will tend 
to suppress such pair field condensation, so we expect the 
supersolid phase with (S) =^ to occur first on increasing 
t z . Such supersolid states have a maximum period of 
2 lattice sites along the principle axes of the triangular 
lattice (as can be seen from the behavior of translations 
in Eqs. QB4[I. so we denote these phases by SS2 (they 
may have doubled or quadrupled unit cells, depending 
upon the orientation of the pseudospin vector). The two 
different ordering patterns for different signs of w 2 are 
shown in Fig. |SJ 




FIG. 8: Possible charge orderings in the SS2 phases. Left: 
S = (100). Right: S = (111). 



two symmetry-unrelated order parameters must change 
in these transitions. We will return to the nature of these 
transitions after first discussing the elementary excita- 
tions of the different phases. 
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FIG. 9: Schematic phase diagram of the Z2 hard spin model. 
Thin lines denote conventional LGW-type transitions, thick 
(red online) lines denote non-LGW quantum critical points. 



C. Elementary excitations 

1. Superfluid phase 

The hard spin model is convenient for describing the 
elementary excitations of the phases discussed above. 
First consider the superfluid. In this case, the elemen- 
tary excitations are simply vortices, and the vortex field 
theory of the previous two sections already gives a de- 
scription of the elementary vortex multiplet, consisting 
in this case of 4 vortex flavors (carrying pseudospin-1/2 
and XY "charge" ±1/2). This should be reproduced by 
the hard spin model. Naively, the "particles" of the hard 
spin model are created separately by the z a and e ±l6 / 2 
fields, so carry only one or the other of pseudospin or XY 
charge. However, in the superfluid region of Fig. |5J the 
Z 2 gauge charges (whose interactions are mediated by 
(Tip) are confined, so the true elementary excitations are 
Z 2 gauge-neutral bound states ~ e ±l9 / 2 z which have 
precisely the appropriate quantum numbers of the z± a 
vortices. 



Putting the different limits of this analysis together 
and making the simplest possible interpolation, we ex- 
pect the phase diagram in Fig. [5] The Mott state may 
be reached from the superfluid in at least three distinct 
ways: by a direct transition described by the continuum 
vortex Lagrangian in the previous section, or via two 
distinct intermediate supersolid phases. The transitions 
from the superfluid to the two supersolids are "conven- 
tional" , i.e. of LGW type, since they are described by 
ordering of the gauge-invariant order parameters ip and 
S. The transitions from the supersolids, by contrast, 
are unconventional. This is clear from the fact that the 
Mott insulator differs from either supersolid by breaking 
more spatial symmetry and by having no off-diagonal 
long range order, i.e. by being non-superfluid. Thus 



2. Mott phase 

Next consider the Mott state. In reality this comprises 
a number of different phases, depending upon the signs of 
the anisotropies w 2l wz- However, near to the superfluid- 
Mott transition, the latter terms are small, and these 
distinct phases are approximately unified into one con- 
tinuous manifold. It is useful to discuss the elementary 
excitations therefore in the same approximation, and af- 
terward describe how they are modified once anisotropy 
is included. From the point of view of the Z 2 hard spin 
model, Eq. H57J) . the Mott state is a Higgs phase, with 
both XY and CP 1 fields condensed. The Z 2 gauge field 
can be regarded, in a choice of gauge, as uniform er^ « 1, 
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in the ground state. If we neglect the possibility of de- 
forming the Zi gauge field in excited states, there are 
then two "obvious" topological excitations - time inde- 
pendent solitons that behave as quantum particles - cor- 
responding to textures in the CP 1 and XY fields. 

First consider the CP 1 field. At spatial infinity, z% a 
must vary slowly in space to maintain minimal action, as 
must Ai^ for the same reason. We may therefore take a 
continuum limit of the first term in Eq. I|57|) . and, taking 
<?iu ~~ * lj one finds the Lagrangian 

C z « ^g s \d^\ 2 + k(C^ - A^) 2 

+^(e^xd„A lX ) 2 (60) 

with g s , k ~ t z , and 

Cn = —^(z*d ll z a -d li z*z a .). (61) 

Minimum action configurations therefore have 9 M S = 
and A^ — C M at infinity. This requires that, at infinity, 
only the phase of the spinor varies, i.e. z a — £ a e l ® , where 
£ CT is a constant normalized spinor, and may depend 
upon the polar angle from the origin. At infinity, then 

C M = d^Q. (62) 

Single-valuedness of z a allows topologically non-trivial 
configurations in which winds by an integer multiple 
of 2ir, whence 

j> dx^Cfj, = j> dx^A^ = 2im s , (63) 

the integer n s being a topological index. Thus the dual 
flux of such configurations is quantized in units of the 
dual flux quantum. Since the dual flux measures phys- 
ical charge, such solitons are particle-like excitations of 
the Mott state with physical integral boson charge n s . 
While at infinity the spinor varies only through 0, (and 
hence the pseudospin S is constant), this cannot hold 
everywhere in space, since such a "vortex" in must 
have a singularity somewhere. If the spinor is assumed 
to vary everywhere slowly in space, so that a uniform 
continuum limit can be taken everywhere, then the sin- 
gularity is avoided by having the associated amplitude 
vanish at the skyrmion's "core" , e.g. in a configuration 
of the form 

Za = f(ry e £ a + y/1 - f(r)va, (64) 

where £ CT and r) a are two normalized orthogonal constant 
spinors, = 0, and f(r) — > 1 as r — > oo, f(r) — ► 

as r — ► (r is the radial coordinate from the skyrmion 
center). The non-collinear variation of z a indicates a non- 
trivial texture of the pseudospin in the skyrmion. Quite 
generally, if z a is analytic, one can show that 

n s = -J- / d 2 r S • d x S x <9„S, (65) 

47T J 



directly relating the skyrmion number to the pseudospin 
texture. Because the pseudospin itself is constant at in- 
finity, it is apparent that the skyrmion has a finite size, in 
the example of Eq. (|64|l determined by the range of sig- 
nificant spatial variation of f(r). This scale is dependent 
upon details of the Hamiltonian within the Mott phase, 
and in general the above description of the spatially- 
varying pseudospin is valid only if this scale is much 
larger than the unit cell of the Mott charge ordering pat- 
tern. Deep in the Mott phase (i.e far from the Mott 
transition), this may not be the case, and in that case 
there is not necessarily any sharp meaning to the pseu- 
dospin texture. In this sense, the skyrmion/ antiskyrmion 
can be considered as adiabatically connected to a simple, 
and patently non-topological, vacancy or interstitial de- 
fect of the Mott "solid" . 

Near to the Mott transition, however, the skyrmion is 
expected to be large, as we now show. The size of the 
skyrmions is determined by the balance of the anisotropy 
energy £ 2 and the interaction energy (e^xd^Ax) 2 ■ A 
rough estimate for the skyrmion size can be obtained by 
a simple dimensional analysis. The anisotropy energy of a 
skyrmion of size A is of the order of u'2|S'| 4 A 2 , where \S\ is 
the unrescaled amplitude of the pseudospin vector order 
parameter. On the other hand, the interaction energy 
is of the order l/e 2 A 2 . The optimal skyrmion size is 
therefore given by: 

A ~ {w 2 \S\ i e 2 )~ 1/i . (66) 

Close enough to the critical point, since |5| becomes 
small, the skyrmion becomes large, and thus develops 
a topological character. 

Now consider the excitations of the XY field. Taking 
(Tip w 1 in Eq. 1)57(1 . the naive topological excitation con- 
sists of winding 9 at infinity by an integer multiple of 
47r - not 2ir, since the cos(A M 6>i/2) is not 27r-periodic. 
These excitations are paired vortices in the 3-sublatticc 
supersolid order parameter ip. They are neutral, since 
there is no coupling to the dual gauge field. Like an or- 
dinary neutral superfluid vortex, they cost a logarithmic 
energy, neglecting the XY anisotropy term W3. When it 
is included, such a (double-strength) ip vortex becomes 
linearly confined, and converts at long distances to an in- 
tersection point of 12 (!) domain walls, the phase winding 
coalescing into these 12 walls radiating outward from the 
"vortex" core. 

Allowing for a non-trivial texture in the Z2 gauge field, 
however, a third kind of excitation is possible in the Mott 
phase. In particular, one may consider a "vison" or Z2 
vortex, around a point around which any line product of 
Ising gauge fields gives —1. This requires the existence of 
a "cut", a ray emanating outward from the vison along 
which Ising gauge fields crossing the ray are taken nega- 
tive (nevertheless, the Z2 flux is non-trivial only through 
one plaquette at the vison core). Such a cut effectively 
introduces anti-periodic boundary conditions for both z a 
and e l6 l 2 across the cut. That such a configuration is 
possible is of course evident from the definition of z a 
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and 9 in Eq. I|47|) . since the apparent discontinuities in 
z a and e ±l9 / 2 do not affect the eiementary z± a fields. 
The anti-periodic boundary condition forces topological 
defects into both the XY and CP 1 fields. In the XY sec- 
tor, it requires the existence of a ±2tt vortex in ip. As 
for the ±47r vortex above, this costs logarithmic energy 
neglecting W3, and degenerates into a linearly confined 
"source" for (in this case 6) radial domain walls. In the 
CP 1 sector, it is slightly less intuitive. One might have 
expected the occurrence of some sort of "half-skyrmion" 
pseudospin texture. However, this is not the case. Anti- 
periodic boundary conditions require the discontinuity to 
persist all the way down from infinity to the vison core. 
There is thus no way to "relax" the winding singularity 
at infinity into a smooth pseudospin texture. Instead, the 
minimal energy configuration has S spatially constant ev- 
erywhere, i.e. has the form z a = , where winds by 
±7r and £ CT is constant everywhere save in some small core 
region of microscopic size. However, the continuum ac- 
tion Eq. I|6U|) still obtains at infinity, so that finite energy 
configurations still satisfy Eq. I|62|) and = C M at infin- 
ity. Hence, these CP 1 "half- vortex" configurations carry 
fractional boson charges ±1/2. So by taking into account 
Z2 gauge vortices, we find a third class of "elementary" 
excitations in the Mott insulator, "half bosons" with a 
texture in the 3-sublattice supersolid order parameter ip. 
These are linearly confined beyond some length at which 
the 6- fold XY anisotropy becomes significant. 

Of the three types of topological excitations discussed, 
it is interesting to note that only the charge ±1 skyrmion 
remains unconfined at the longest scales. 



3. SS3 phase 

Let us now turn to the SS3 phase, which is described by 
Eq. H57|) at large tg. As discussed above, the ground state 
in this limit can be regarded as a Z2 Higgs phase, with 
the z a field uncondensed. Hence there are two types of 
topological defects: "paired" XY vortices in ip, in which 
9 winds by a multiple of An and single XY vortices, in 
which ip winds by ±27r, accompanied by a vison. Both 
cost logarithmic energy at short scales, crossing over to 
linear confinement as do similar excitations in the Mott 
state. Finally, there are the CP 1 "particles" created by 
the z a field, which can be considered to propagate coher- 
ently since the Z2 gauge field is in a Higgs phase. The 
single XY vortices have a statistical interaction with the 
CP 1 quanta, but this does not lead to significant effects 
upon the CP 1 particles since the XY vortices are anyway 
linearly confined. The CP 1 quanta still carry unit dual 
gauge charge, and so should be regarded as the physical 
superfluid vortex excitations of the supersolid. 



4. SS2 Phase 

The pseudospin vector order parameter S = z*tz is 
condensed in the SS2 phase. However, it is not a Higgs 
phase for the z a fields, since, for instance, it is still super- 
fluid, i.e. the dual gauge field remains gapless. Thus in 
the SS2 phase Z<i quanta are strongly confined, and the 
elementary excitations must be Z2 singlets. One class of 
excitations are skyrmions in S. They do not carry any 
well-defined charge since the boson number conservation 
symmetry is anyway broken in the supersolid. The other 
quanta are physical vortices, which are bound states of 
z a and e ±l6 / 2 particles, essentially the original z± a vor- 
tices of the superfluid. However, because of the broken 
pseudospin symmetry in this phase, there is a preferred 
pseudospin polarization, and the two a components of 
the vortex spinor (choosing a quantization axis along S) 
are no longer energetically equivalent. Therefore there is 
only a two-fold rather than four-fold low-energy vortex 
multiplet £± ~ z±q, taking a = as the lower-energy 
spinor. 



D. Supersolid-Mott quantum critical points 

1. SS3-Mott transition 

The supersolid to Mott insulator QCPs are manifestly 
not of LGW type, if they occur at all as continuous 
transitions. Nevertheless, at least some aspects of these 
QCPs can be straightforwardly analyzed by application 
of Eq. J37J). Consider first the SS3-Mott transition. It 
is useful to approach the transition first from the SS3 
phase. Since it can be regarded as a Z2 Higgs state, 
this is particularly simple. In particular, we can treat 
(7^ « 1 as a constant at low energies. The e ±l6i ^ 2 fields 
can be regarded as condensed, and moreover at low en- 
ergies there are no associated gapless Goldstone modes 
due to the 6-fold anisotropy term w 3 in Eq. (|4*1*1) . Hence 
at low energies, the only important fluctuations are those 
of the CP 1 spinor and the dual gauge field. The natural 
critical theory, kept lattice regularized for simplicity, and 
neglecting for the moment pseudospin anisotropy, is thus 
just 

C-nccpi = —t^i^c lA ^z* ic! Zi + ^ a + — -tj (efj,v\A u Ai\) 2 . 

(67) 

This is the Non-Compact CP 1 (NCCP 1 ) theory studied 
numerically in Ref. Il5l and believed to represent a dis- 
tinct universality class of critical phenomena. It also has 
a more intuitive interpretation: it describes the behavior 
of the 2 + 1 = 3 dimensional 0(3) transition associated 
with S when "hedgehog" defects in S are completely sup- 
pressed in the partition function. Let us now see how 
this is understood approaching the QCP from the Mott 
side. In this phase, the pseudospin vector is ordered, 
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but fluctuates more and more strongly as the SS3 phase 
is approached. One would expect skyrmion defects to 
become more and more prevalent as fluctuations in the 
pseudospin increase. As described above in Sec. IIV C 21 
however, skyrmions in the Mott state carry physical bo- 
son charge, so boson number conservation requires that 
skyrmion- number (n s in Eq. H65|) ~) must also be con- 
served. Happily, skyrmion number-changing events are 
exactly the hedgehog defects of the 0(3) model, so we 
see that charge conservation causes this transition to be 
of the NCCP 1 type. 

Eq. I|67|l and the subsequent discussion neglect pseu- 
dospin anisotropy. While it is quite likely such anisotropy 
terms are irrelevant the NCCP 1 fixed point, this requires 
further study. Using the hard-spin PSG transformations, 
Eqs. (|C2|) . the leading anisotropy terms can be shown to 
be 

Z'nccpi = W2 J2(S a ) 4 + w' 2 (lm^)S x SyS z .(68) 

a—x,y,z 

The latter term is allowed by symmetry, but vanishes for 
W3 > 0, in which case tp 3 is purely real. Note that the 
perturbations w 2 , w' 2 in Eq. l|6^|) are 8 th and 6 th order in 
the CP 1 fields, respectively, so it is quite plausible that 
both are irrelevant at the NCCP 1 point, though clearly 
this is most likely for W3 > 0, when the w' 2 term is absent. 

A further complication in the case W3 < is that 
the non- vanishing -0 order parameter in this case breaks 
particle-hole symmetry C (actually the supersolid with 
W3 > also breaks C, but preserves the combination 
C o T 2 o i?2ir/3 Idi Id 2 i which is sufficient). Since a 
supersolid, like a superfluid, is compressible, this has the 
difficulty that it implies a non- vanishing deviation of the 
spatially averaged density from half-filling (working at 
fixed chemical potential chosen to maintain particle-hole 
symmetry - i.e. zero Zeeman field in the XXZ model). 
In the canonical ensemble, fixing the average density at 
/ = 1/2, it implies phase separation. Formally, this is de- 
scribed in the dual theory by the allowed coupling term 

C" = \(lmip 3 )(A x A y -A y A x ), (69) 

since the physical density is the dual magnetic flux. This 
indeed leads to a density deviation from half-filling away 
from the NCCP 1 critical point, since it leads to a min- 
imum of Cnccpi + £" with non-zero Sn ~ (A x A y — 
A y A x )/2ir. As the Mott transition is approached, how- 
ever, the system becomes increasingly less compressible, 
and the compressibility certainly vanishes when super- 
fluidity does. It is not entirely clear to us how this is 
resolved - the complications are similar to (but more 
difficult due to the pseudospin structure) those occur- 
ing in the theory of the normal-superconducting thermal 
phase transition in a three-dimensional superconductor 
in a weak applied external field Hml It is possible that, at 
fixed chemical potential, the NCCP 1 critical fixed point 
is "weakly avoided" at long scales by this effect, most 



likely by introducing a narrow region of an "SS6" phase 
- a supersolid with the same symmetry as the Mott in- 
sulator but with ODLRO - between the Mott insulator 
and ferrimagnetic SS3 state. Working at fixed density 
/ = 1/2, one expects to pass through the NCCP 1 point, 
which coincides with the critical endpoint of the phase 
separation region. 

2. SS2 to Mott transition 

As indicated in Sec. IIV C 41 the SS2 phase should be 
thought of as a state in which (S) 7^ 0, but vortices 
themselves are not condensed. It is not, however, a 
Higgs phase of Eq. I|57|l . Moreover, the important ele- 
mentary excitations of this phase are just vortices, which 
are bound states of the hard-spin fields. Therefore it is 
advantageous to return to the original soft-spin vortex 
formulation, and procede by just adding the term 

4j S2 = -A<S>-(S+ + S_), (70) 

where of course the S± fields on the right should be 
understood as composites of z± a . As discussed in 
Sec. IIV C 41 this splits the 4-fold vortex multiplet into 
two 2-fold multiplets. Taking (S) along (100) (we will 
not discuss the (111) case in any detail, but it is simi- 
lar), we obtain the scalar low energy fields z±, defined 
by z± a = z±rj a , with t x tj — +rj (other orientations are 
solved by the obvious generalization). By considering the 
residual symmetry operations of this SS2 state (see Ap- 
pendix Eli we may thereby derive the continuum action 
for these two fields: 

£= E [\^~^)z a \ 2 + S \z a \ 2 ] +u(z* a z a f (71) 

a=± 

-v\z + \ 2 \z_\ 2 + \lm(z* + z_) 6 + -L^^Aa) 2 . 

Here we have kept the 12 th order A term because it is 
the lowest order term which breaks the staggered U(l) 
symmetry. 

Remarkably, Eq. (|71|l is extremely similar to the 
NCCP 1 Lagrangian, differing mainly in that the SU(2) 
symmetry of that theory is here reduced by the v term 
to U(l). For v > 0, corresponding to the easy-plane case, 
it is the continuum theory for the "deconfined quantum 
critical point" of Refs. that describes the superfluid to 
VBS transition on the square lattice, with the modifica- 
tion that the "clock" anisotropy A is here 6-fold rather 
than 4-fold. This theory, neglecting the irrelevant A term, 
is self-dual at the critical point, and can alternatively be 
formulated as a theory of the bose condensation of charge 
±1/2 fractional bosons. These are just the half-boson ex- 
citations discussed in Sec. IIV (j2l on the Mott state, which 
carry a direct U(l) gauge charge. 

On examining the charge ordering pattern in Fig. [3] an 
interesting question arises. The symmetry of the (100) 
SS2 state is already consistent with a very simple half- 
filled Mott insulator, consisting of stripes of charge on 
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alternate lines of sites (along principle axes of the trian- 
gular lattice). So it is perfectly conceivable that in some 
models, one could have a transition from the SS2 super- 
solid to a Mott insulator with the same symmetry as the 
SS2 phase. One would expect this to be an XY transi- 
tion, since only the superfluid U(l) symmetry is broken 
across the transition. Why does our theory not predict 
this simpler scenario? 

Firstly, we note that deconfined quantum criticality 
for the SS2 to Mott insulator transition studied above 
is perfectly consistent, since the Mott insulator in ques- 
tion is not the one with the same symmetry as the SS2 
phase. The question remains why we do not see that 
possibility as well. Our interpretation is that, by start- 
ing with the dual field theory for the z± a vortices, we 
have chosen a restricted set of vortex modes (this partic- 
ular multiplet), which describes the natural instabilities 
of an isotropic triangular lattice superfluid. Though we 
have lowered the symmetry already in the SS2 phase, we 
have presumed the low energy excitations in this phase 
should be taken same vortex multiplet. If the spatial 
symmetry breaking present in the SS2 phase were taken 
strong, this might not be a good assumption, and states 
origination from other vortex multiplets could cross the 
z± states in energy, and lead to instabilities to different 
Mott states and also different critical behavior. We leave 
an exploration of this idea for future work. 



3. SS2 to spin liquid transition 

On passing from the SS2 phase to the Z2 spin liquid in 
Fig. a vortex pair field z a z a i must condense. Because 
in the SS2 phase, the pseudospin rotational symmetry 
is already broken, we expect the pair field composed of 
two spinors aligned along the S axis to describe the con- 
densate. This is a one-component field with dual U(l) 
gauge charge 2. Hence we expect this transition to be 
described by a massless charge scalar coupled to a non- 
compact U(l) gauge field. This is just dual to the XY 
model, so this can be viewed as an XY transition. In 
more physical terms, on passing from the Z2 phase to 
the SS2 phase, the half-boson excitations of the spin liq- 
uid condense. It is clear from this description that the 
symmetry of the spin liquid is the same as that of the 
SS2 phase. 



4- Spin liquid to Mott transition 

In the Z2 spin liquid, both S and the vortex pair field 
are condensed. It has, however, the symmetry of the SS2 
phase. It can be viewed as the Higgs phase of the 
gauge field. To pass to the Mott state, which does not 
have topological order, vison excitations must condense. 
The e ±%e l 2 particles play the role of the vison. This can 
be seen, e.g. from the fact that these particles have a sta- 
tistical interaction with the half-bosons in the Z2 phase, 



which are 7r-flux tubes in Ai^ . The non-standard feature 
is that the visons also carry space group quantum num- 
bers - since the e ±%B l 2 operator is a "square root" of the 
SS3 order parameter (the 9 transformations are given in 
Appendix ICT) . Actually this transition can be understood 
from Eq. (|57J) simply by "freezing" ~ 1 and treating 
the Zfj fields as condensed. It is thus clearlyan XY tran- 
sition, with either 6-fold or 12-fold "clock" anisotropy (it 
is doubled since "half" the SS3 order parameter is con- 
densing), for S along (100) and (111), respectively. 



E. Anti-parallel pseudospins 

Here we briefly sketch the results of an analogous study 
of the antiparallel pseudospin case. Using the parame- 
terization in Eq. I|48[> to define the hard-spin degrees of 
freedom, and gauging the Z2 redundancy, the appropri- 
ate hard spin model is 

£>zl = -tzVinZio-Zi+fKr ~ cos(A M 0i/2 - A ifl ) 

+^{e lxu \KA i \f. (72) 

Note that, in this case, the dual gauge field is coupled to 
the XY and not the CP 1 de gree of freedom. 
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FIG. 10: Schematic phase diagram for antiparallel pseudospin 
vectors. Thick and thin lines have the same meaning as in 
FigEl 

Study of various limits and interpolation determines 
the phase structure of Eq. 1|72[) . which is schematically 
shown in Fig.^31 In addition to the superfluid and Mott 
states occuring when both t z ,tg are small and large, re- 
spectively, there is a supersolid phase, SS12 (it has a min- 
imal unit cell of 12 sites), and a Z2 spin liquid insulating 
phase. Unlike in the parallel pseudospin model, the SS12 
supersolid has the same symmetry as the nearby Mott 
insulator. This is a direct consequence of the fact that 
the dual gauge charge is in this case carried by the U(l) 
degree of freedom e 10 , and this does not carry any space 
group quantum numbers. For the same reason, the Z2 
spin liquid phase at large t$ does not break any symme- 
tries. Note that of the four phases in Fig.^J only the Z2 
spin liquid is "exotic" , i.e. has an underlying topologi- 
cal order and unconventional excitations not captured by 
any local mean-field theory and order parameters. 
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The excitations are as follows. In the Mott state there 
are neutral skyrmions (since the z a fields carry no dual 
gauge charge), confined charge ±1/2 excitations (visons 
bound to half-vortices in z a ,9), and charge ±1 excita- 
tions that can be viewed as ±4-7r vortices in 6. Because 
9 is coupled to Ai^, the latter cost finite energy, and are 
clearly adiabatically connected to vacancy/interstitials in 
the Mott state. The SS12 phase has the same skyrmion 
textures, but no well-defined charge excitations since it 
is a superfluid. Instead, it has a single physical vor- 
tex/antivortex excitation created by e ±l6 / 2 . The "trivi- 
ality" of the vortex multiplct is consistent with the bro- 
ken symmetry of the supersolid, whose enlarged unit cell 
contains on average an integer number (6 in the sim- 
plest case) bosons. The Z2 spin liquid has physical bo- 
son charge ±1/2 excitations {2tt vortices in accompa- 
nied by a "vison" in ct^) and physical "vison" excitations 
(created by e ±l8 ^ 2 and z a ) which carry spatial quantum 
numbers. 

The direct transition from superfluid to antiparallel 
Mott state is described by the continuum action of the 
previous section. The superfluid-SS12 and the SS12-Mott 
critical points are also conventional, since each is char- 
acterized by a change in a single order parameter. The 
superfluid-SS12 transition is described by an LGW the- 
ory for the d-vector, while the SS12-Mott transition is 
simply an XY transition for the superfluid order pa- 
rameter. The superfluid- Z2 transition is also an XY 
transition, which can be understood as a condensation 
of charge ±1/2 "half-bosons" (in principle this changes 
universal amplitudes from the conventional superfluid to 
integer- filling Mott transition, which is also XY-like^). 
The Z2 to Mott transition is modeled by the CP 1 action 
with no gauge field, which has the physical interpretation 
of modeling vison condensation. We have not attempted 
to consider the effects of various anisotropies on these 
transitions. 



V. DISCUSSION 

In this paper we have presented a phenomenological 
dual vortex theory of the interplay between Mott local- 
ization and geometrical frustration for interacting bosons 
at half-filling on the triangular lattice. This approach re- 
veals a variety of novel quantum phases and phase tran- 
sitions which may occur if the superfluid and Mott in- 
sulating states occur in close proximity to one another 
in phase space. Of particular interest are the continuous 
superfluid-Mott insulator transition predicted by mean 
field theory, the two supersolid phases, and the occur- 
rence of the recently-discovered NCCP 1 critical universi- 
tality class at the 3-sublattice supersolid to Mott insula- 
tor transition. In this discussion, we will provide a more 
direct physical picture of some of these phenomena, and 
address the prospects of observing them in simple micro- 
scopic boson or spin models. 

A useful starting point for the discussion is the recent 



demonstration that a supersolid phase indeed occurs in 
the simplest spin-1/2 XXZ model, 

Hxxz = ^ —J±(S^Sj + SfSj) + J z S*Sj, (73) 

with ferromagnetic XY and antiferromagnetic Ising ex- 
changes (Jj_, J z > 0) (equivalently, hard-core bosons with 
nearest neighbor repulsion) on nearest-neighbor links of 
the triangular lattice ^LiLiS This model was shown to 
be in a 3-sublattice SS3-type phase for J z > 5J±, and 
this phase persists up to and including J z — 00. A 
number of features of the numerical results on the su- 
persolid at large J z are notable. First, although super- 
fluidity survives, it is extremely weak, as characterized 
by the superfluid stiffness, which is approximately 250 
times smaller in the large J z supersolid than in the pure 
XY model ( J z — 0). Second, it is exceedingly difficult to 
distinguish numerically on even relatively large lattices 
between the two different types of SS3 charge ordering 
patterns. Ref. llCll was unable to distinguish them numer- 
ically by direct measurement of boson density correlation 
functions, while Ref. ^ claimed to do so, but only for 
large lattices of 18 x 18 sites with a very small signal. 
Furthermore, a deviation of the density from half-filling 
is expected 10 in the "ferrimagnetic" SS3 phase identi- 
fied in Ref. ^] (which corresponds in our theory to the 
one with W3 < 0), but appears to be exceedingly minute 
if observable at all computationally. Apparently there 
is very little splitting energetically between the two SS3 
states, even at J z — 00, a point at which there is no in- 
trinsic small parameter in the microscopic Hamiltonian 
- the effective Hamiltonian is simply the XY exchange 
projected into the Hilbert space spanned by the mani- 
fold of classical Ising antiferromagnetic ground states on 
the triangular lattice. 

The present theory offers a partial explanation for 
these puzzling observations. We interpret the weakness 
of superfluidity as evidence that the system is in close 
proximity to a Mott insulating state. If so, our dual 
vortex field theory, which is built around a superfluid to 
Mott insulating transition, should apply. The tiny energy 
splitting between the two SS3 states is then understand- 
able: the term which dictates this splitting, W3 is 12 th 
order in the basic z± a vortex fields, and clearly strongly 
irrelevant at the superfluid-Mott QCP. Furthermore, the 
smallness of any spontaneous density deviation from 1/2- 
filling, which as indicated above is expected for the "fer- 
rimagnetic" SS3 state, is also expected, since the density 
deviation must vanish as the Mott state, which has den- 
sity of exactly 1/2 and is incompressible, is approached. 
This deviation also vanishes at the transition from this 
state to the superfluid (as seen in LGW theorjiiS), so it 
is likely small throughout the ferrimagnetic phase. 

Given these arguments for proximity to the Mott state, 
it seems likely that only a small perturbation of the 
XXZ Hamiltonian may be required to push it into a 
Mott phase, and in so doing observe the very interest- 
ing NCCP 1 criticality at the SS3-Mott quantum critical 
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point. Let us try to develop a more physical picture of 
this transition. We will begin by providing a cartoon un- 
derstanding of th e SS3 phases. In the case W3 < 0, as 
discussed in Refs. I 111 1 2t the ferrimagnetic SS3 state can 
be understood crudely by first forming a "solid" of bosons 
with a density of 1/3, with one boson occupying each of 
the sites of one of the three \/3 x y3 triangular sublat- 
tices of the original lattice. This leaves, at / = 1/2, a 
density of 1/4 boson per the remaining sites, which form 
a honeycomb lattice, with the "solid" bosons in the cen- 
ters of the honeycomb plaquettes. The SS3 state can 
be viewed as a superfluid of these remaining bosons (al- 
ternatively, one can make the same construction with 
holes replacing bosons, leading to a different but equiv- 
alent state - which illustrates the spontaneously broken 
particle- hole symmetry of the ferrimagnet). In the op- 
posite case W3 > 0, there is no spontaneous density po- 
larization, and the "sublattice magnetizations" take the 
values (Sf) — (m, — m,0) on the three inequivalent sites. 
We will call this the "antiferromagnetic" state because of 
the exactly opposite Ising moments on two of the three 
sublattices. This state can be understood by a similar 
cartoon. In particular, take a honeycomb sublattice of 
the triangular lattice, and on this sublattice form a 1/2- 
filled Mott insulator of alternating empty and occupied 
sites (an Ising Neel state in spin language). The sites at 
the centers of the honeycomb plaquettes form a V3 x \/3 
triangular sublattice, and we put the remaining bosons 
into a half-filled superfluid on this sublattice. These car- 
toon pictures would clearly tend to favor the ferrimag- 
netic state in the XXZ model, since bosons must hop 
(presumably virtually) between second neighbor sites to 
stabilize the triangular superfluid in the antiferromag- 
netic state. This immediately suggests that, to study 
the antiferromagnetic supersolid, one needs only to add 
a second neighbor XY exchange to the XXZ model, 

H' = -J' ± J2 S i S I + S ^ S J- ( 74 ) 

With these pictures of the SS3 phases in hand, it is 
natural to view the transition to the Mott insulator as a 
"crystallization" of the superfluid sublattices of the su- 
persolid. Indeed, we find that this provides a simple 
physical picture consistent with the symmetries of the 
appropriate Mott phases. Consider first the transition 
from the ferrimagnet. Here we have a 1/4- filled honey- 
comb lattice of bosons, which is to undergo a superfluid 
to Mott insulator transition. At 1/4-filling, these bosons 
are unlikely to form a simple "crystalline" Mott insula- 
tor, since they would be too widely separated to substan- 
tially interact. Instead, more likely Mott states are va- 
lence bond solids, in which the bosons resonate between 
two or more sites (within still-localized wavefunctions). 
The most natural candidate is a "columnar" valence bond 
solid (VBS) state, in which alternating columns of bonds 
are occupied by one boson. Indeed, the state predicted 
by the vortex mean field theory for W3 < and W2 < 
has exactly the symmetries of the columnar valence bond 



solid, see Fig. ^] One may also convince oneself of the 
validity of the columnar VBS picture by counting the 
number of distinct Mott states. Fix the location of the 
V3 x \/3 superlattice and hence the honeycomb sublat- 
tice. One can then place the valence bonds along columns 
parallel to any of the three principle axes (of the honey- 
comb, which has principle axes halfway between those of 
the underlying triangular lattice), and for each such ori- 
entation, they may lie on even or odd columns. Hence one 
expects six states. This is precisely the number of dis- 
tinct choices of vector pseudospin along the ±x, ±y, ±z 
axes. 




FIG. 11: Density plot for the ferrimagnetic solid (left) and 
the corresponding cartoon (right). In the density plot, the 
honeycomb sublattice is clearly visible. In the cartoon, the 
honeycomb lattice on which the superfluid lives is overlaid in 
green, with the frozen bosons in blue, and the valence bonds 
that form on passing from the ferrimagnetic SS3 phase to the 
ferrimagnetic solid are drawn as red ovals (color online). 

Next consider the antiferromagnetic supersolid. We 
have an effective \/3 x \/3 triangular sublattice of bosons 
at half-filling, living in the plaquette centers of an "an- 
tiferromagnetic" bose solid on the honeycomb lattice. In 
the absence of this surrounding solid, the triangular sub- 
lattice would have, scaled up to its size, all the same 
symmetries as the original triangular lattice. The for- 
mation of a Mott insulator on this triangular sublattice 
would thus naively appear to be just as formidable a 
problem as the original one. The staggered solid, how- 
ever, breaks Id 2 and C, preserving only the combination, 
Id 2 o C. This means, were we to repeat the PSG anal- 
ysis for a dual vortex theory of the -\/3 x y/3 sublattice 
bosons, there is no symmetry to prevent the Ising order 
parameter $ = |z+| 2 — \z~\ 2 from appearing as a term 
in the action, breaking the a = ± "flavor" degeneracy 
of the vortex multiplet. Clearly in such a theory, then 
only one of the two flavors will condense, and only those 
phases in which one pseudospin is non-zero will appear. 
Incidentally, one may identify therefore the z a spinor in 
the hard-spin action with the lower-energy vortex flavor 
in this sublattice theory. The phases with only a sin- 
gle pseudospin condensed are "staggered" phases of the 
original bosons. The simplest of these is just a "columnar 
crystal" , in which bosons on alternating columns of trian- 
gular lattice sites (along some principle axis) are occupied 
and unoccupied. Superimposing this upon the surround- 
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ing antiferromagnetic honeycomb lattice, one remarkably 
obtains an ordering pattern again identical in symmetry 
to the original mean-field solid with W3 > and W2 < 0! 
The counting of states is also the same as for the ferri- 
magnetic case above, since the "columns" have the same 
set of orientations and have only shifted from bonds to 
sites, once again in agreement with the configurations of 
S. 

One can easily extend these constructions to the cases 
with pseudospin along (111). For brevity, we relegate 
this to Appendix El 




FIG. 12: Density plot for the antiferromagnetic solid (left) 
and the corresponding cartoon (right). In the cartoon, the 
triangular lattice on which the superfluid lives is overlaid 
in green, with the crystallized bosons/vacancies forming the 
Mott insulating honeycomb lattice in blue (filled and empty 
circles). The additional bosons/vacancies that solidify on 
passing from the antiferromagnetic SS3 phase to the antifer- 
romagnetic solid are drawn as red squares (filled and empty, 
respectively) . 

The cartoons of the ferrimagnetic and antiferromag- 
netic supersolids/solids not only reproduce the symme- 
tries of the phases, but also can be directly used to ob- 
tain the critical NCCP 1 theory of the SS3-Mott transi- 
tion. In the ferrimagnetic case, the cartoon picture is 
a superfluid-Mott transition on a 1/4- filled honeycomb 
lattice of bosons. In Appendix IF1 we sketch how a dual 
vortex analysis of that "cartoon" problem directly recov- 
ers the appropriate NCCP 1 theory for this case. In the 
antiferromagnetic similar argument has already 

been sketched for the cartoon of an half-filled triangular 
V3 x V3 sublattice in the background of the honeycomb 
antiferromagnetic Mott state, which removes one of the 
two pseudospinors from the theory, likewise recovering 
the NCCP 1 model and the appropriate anisotropy. 

It thus seems of considerable interest to pursue small 
perturbations of the XXZ model which might take it 
into a Mott phase, and between ferrimagnetic and an- 
tiferromagnetic states. To make the latter transition, 
from the ferrimagnetic to antiferromagnetic SS3 super- 
solid, we can make a simple recommendation: the inclu- 



sion of weak second-neighbor ferromagnetic exchange, as 
in Eq. I|74|l . should favor the antiferromagnetic state by 
allowing bosons to better move on the triangular sublat- 
tice. The antiferromagnetic states seem slightly more fa- 
vorable than the ferrimagnetic ones theoretically as can- 
didates for observing the NCCP 1 criticality, because the 
latter preserve particle-hole symmetry in the supersolid, 
and have weaker pseudospin anisotropy terms. Under- 
standing what terms microscopically would be needed to 
push the XXZ model toward the insulators is a nontriv- 
ial problem. However, it seems clear that introducing 
further-neighbor Ising exchange will lead to some Mott 
insulating states, when J z is large, since they spoil the 
degeneracy of the classical Ising manifold. Exploration 
of this problem numerically is tempting. 

Analytically, it is quite interesting how naturally the 
dual vortex approach leads to the observed supersolid 
states. Coming from the Mott insulator with paral- 
lel pseudospins, the deconfined charge excitations are 
skyrmions. If one imagines lowering the charge gap in 
the insulator, it is most natural that these should con- 
dense, destroying the pseudospin order of the insulator 
and simultaneously initiating superfluidity. While other 
scenarios such as the direct Mott insulator to superfluid 
transition are possible theoretically (and described ana- 
lytically in this paper), they entail a non-trivial screening 
so that confinement of other (fractionally charged) exci- 
tations of the Mott state becomes progressively weaker 
as the transition is approached. Thus the SS3 supersolid 
is probably the most intuitive "neighbor" of the Mott 
state, and with this understanding its occurence in the 
strong interaction limit of the XXZ model is no longer 
surprising. It is interesting to speculate as to whether 
it might be possible to formulate the boson-vortex du- 
ality used heavily in this paper directly in this strong- 
coupling region, i.e. in the classical Ising antiferromag- 
netic Hilbert space. Such a formulation could potentially 
provide a quantitative method for attacking problems of 
Mott transitions with strong frustration. 
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The PSG transformations in momentum space are given by: 
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Using these general forms for the PSG in momentum 
space, it is straightforward to obtain Eqs. 1)12113(1 of the 
main text. The remaining PSG transformations for the 
multiplets are: 



It is possible that 771(a) = 772(a) holds at a general 
filling, not just at / = 1/2, but it is not obvious. 

For the specific case of / = 1/2 (p = 1, q = 2), we can 
impose invariance under the particle-hole transformation, 
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Imposing invariance under the rotation and reflection 
operations above upon the quartic action C\ taken in the 
form of Eq. (|14fl yields the set of conditions 
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for g odd, and 



when g is an even integer. Appying ^2^/3 to one 
finds that 
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Applying , one finds that 773(a) = 0. The remaining 
phase factor ?y 2 (a, /) is difficult to find analytically in the 
general case. In the case / = 1/2 considered in detail in 
the text, it is, however, equal to 771(a), so collecting these 
results, 



1 9-1 

v a/3 _ _|_ \ ^ 



n ,afi , ,mn' —nm' 



Im'n' ~ / 1 i' 
m,n— 



m,n— 
9-1 



771(a) = 7y 2 (a) = -Tra/12, 773(a) = 



for / = 1/2. 

(A5) 



for q even. 
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APPENDIX B: PSEUDOSPIN 
TRANSFORMATIONS AT / = 1/2 

For convenience, we give the transformation proper- 
ties of the spinor z aa vortex helds. With the pseudospin 
index suppressed, we find 
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(1,0,-1)/V2. 



Next we tabulate the transformation properties 
pseudospin vector order parameters: 
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APPENDIX C: HARD SPIN TRANFORMATIONS 

Here we give the transformation rules for the hard-spin 
theory with parallel pscudospins. Defining 



z ± = ze^l\ 



(CI) 



we obtain 
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(C2) 

The SU(2) matrices u are defined in Eqs. i|B2() . We have 
used the freedom to redefine any of these operations by 
a global U{\) gauge transformation, which in this hard- 
spin limit corresponds to a phase rotation of z. 



APPENDIX D: RESIDUAL SYMMETRIES OF 
THE SS2 SUPERSOLID 

The residual symmetry operations of this SS2 state 
(which are not broken by the pseudospin vector order) are 
generated by Tx , R 27T / 3 ° Idi , C o I d2 o I dl , and T 2 oI d2 oI dl . 
Their action on the z± vortex fields can be obtained from 
the definition z± a = z±r\ a and Eqs. i|Bl|) : 



Tx : z± -> ±e^ 9 z±, 



R 2 Tr/3 ° hi 



Z- 



C o I d2 o I dl : z± -> < 
T 2 o I d2 o I dl : z± - 

These lead directly to Eq. (J2TJ 
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APPENDIX E: TOY MODEL WAVEFUNCTIONS 
WITH (111) PSEUDOSPIN 

Here we provide toy model wavefunctions for the Mott 
states with w 2 > and parallel pseudospins, so that 
S + = S_ lies along say the (111) axis. First consider 
W3 > 0, in which we wish to construct a state appro- 
priate to a 1/4-filled honeycomb lattice. In this case, 
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the boson Green's function. 



FIG. 13: Density plot for the ferrimagnetic solid (left) and the 
corresponding cartoon (right) with (111) pseudospin ordering. 



a wavefunction with the correct symmetry is shown in 
FigED I n this figure the bosons denoted by squares are 
on the honeycomb sublattice of the triangular lattice, and 
become superfluid in the SS3 phase. 

Now consider W3 < 0, in which case we must consider 
the half-filled triangular lattice embedded in the antifer- 
romagnetic honeycomb Mott insulator. The density plot 
for this state is shown in Fig ll4l Focusing on the trian- 
gular lattice sites, we note that the amplitude on these 
sites has the same 2x2 periodicity as the W3 > state 
above, with the sites on this sublattice having one am- 
plitude, and the remaining sites another. However, it is 
inconsistent with half-filling to simply put one boson on 
the sublattice sites, or on the other three sites. Moreover, 
the state is rather constrained by the existence of a num- 
ber of centers of three- fold rotations and reflections. The 
simplest Mott wavefunction we constructed that satisfies 
all symmetry properties is actually not quite a product 
state of the form of Eq. H25[l . It is, however, clearly an in- 
sulating state. To construct it, assume the 2x2 sublattice 
sites are empty. Then the remaining sites form a kagome 
lattice. The kagome lattice is composed of corner-sharing 
triangles, or two orientations (pointing "left" and "right" 
in the figure). We act with creation operators to place 
one boson on each triangle in a uniform superposition. 



*=IMl+kA2+^ 3 )|0> 



(El) 



A 



Since the triangles overlap, there will actually be ampli- 
tude to find more than one boson per triangle. On av- 
erage, the number of bosons per unit cell of the kagome 
lattice is then 2, since there is one up and one down tri- 
angle per unit cell. There are four triangular lattice sites 
for each such unit cell (3 kagome and one central empty 
site), so this is a state at 1/2-filling, which manifestly has 
all the symmetries of the Mott state in Fig^J Moreover, 
it is an insulating state as required, since one can readily 
show there are only very short-range correlations, e.g. in 




FIG. 14: Density plot for the antiferromagnetic solid (left) 
and the corresponding cartoon (right) with (111) pseudospin 
ordering. 



APPENDIX F: HONEYCOMB LATTICE AT 
1/4-FILLING 

The vortex PSG's on the honeycomb lattice were 
worked out in Appendix B of the first paper in Ref. 0. 
They are sufficient to develop a dual vortex theory of 
the superfluid-Mott transition on this lattice. For the 
case of / = 1/4, there are two vortex flavors, ipo,ipi, 
which can be combined into a spinor ip a . We consider 
the restrictions upon the vortex Lagrangian by transla- 
tions and rotations. Suppressing the indices and using 
Pauli matrices T aa > to span the spinor space, we have: 
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der these PSG operations, one may readily write down a 
continuum Lagrangian in terms of p: 

Aioney = ^ [\ (<9 M ~ iA^)<p a | 2 + s\p a | 2 ] + u(p* a Lp a ) 2 



<T=0,1 



+ -^{^^u\d v A\) 2 + C a , 



(F2) 



where SU(2) symmetry of <p a is first broken by the term 



C a = \S X S V S Z 



(F3) 



with S = Lp* a T aa ip> a . Eq. (|F2I) indeed has the global 
SU(2) and U(l) gauge symmetry of the NCCP 1 model, 
and C a recovers the leading anisotropy term in Eq. I|68() 
for the ferrimagnetic case. 
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